Production and polarization of the A^ and the 
Q^ ■ charm of the proton * 

;-i ! J. C. Anjos G. Herrera ^ J. Magnin ^ 

^ ■ F.R.A. Simao 

'!;;j- ' Centre Brasileiro de Pesquisas Fisicas 

I Rua Dr. Xavier Sigaud 150 

^ ! 22290-180, Rio de Janeiro, Brazil 

^: 

vn ■ 

(N : PACS 13.60.Rj (Baryon production) 

j^ I 13.87.Fh (Fragmentation into hadrons) 

0> ' 14.20.Lq (Charmed baryons) 

PhI Abstract 

rr^' We propose a two-component model involving the parton fusion 

r~| , mechanism and recombination of a ud valence diquark with a sea c- 

^ I quark of the incident proton to describe A+ inclusive production in 

pp collisions. We also study the polarization of the produced A+ in 

the framework of the Thomas Precession Model for polarization. We 

Cd I show that a measurement of the A^, polarization is a sensitive test of 

its production mechanism. In particular the intrinsic charm model 
predicts a positive polarization for the Ac within the framework of the 
Thomas Precession Model, while according to the model presented 
here the Ac polarization should be negative. The measurement of 
the Ac polarization provides a close examination of intrinsic charm 
Fock states in the proton and gives interesting information about the 
hadroproduction of charm. 
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1 Introduction 

The production mechanism of hadrons containing heavy quarks is not well 
understood. Although the fusion reactions gg -^ QQ and qq —>■ QQ are sup- 
posed to be the dominant processes, they fail to explain important features in 
hadro-production like the leading particle effects observed in D^ produced in 
71^ p collisions [jl|, A+ produced in pp interactions 0, ||, Q and others baryons 
containing heavy quarks [0, the J/\E' cross section at large xp observed in 
Tcp collisions 0, etc. 

The above mentioned effects have been explained using a two-component 
model [0 which consists of the parton fusion mechanism, calculable in per- 
turbative QCD, plus the coalescence of intrinsic charm ^j. 

Here we present an alternative mechanism, namely the conventional re- 
combination of valence spectator quarks with a c-quark present in the sea of 
the initial hadron. 

We show that this mechanism can explain the production enhancement at 
large Xp observed in charmed hadron production in hadronic interactions. 

We describe the A^ production in pp interactions with a two-component 
model consisting of the recombination of a diquark ud with a c-quark from 
the sea of the incident proton plus the usual parton fusion and fragmenta- 
tion mechanism. We compare the obtained results with those of the intrinsic 
charm two-component model and with experimental data. 

We also study the polarization of A+ using the Thomas Precession Model 
(TPM) for both the intrinsic charm model and the conventional recombi- 
nation two component model presented here. We show that the polarization 
prediction of this two models is totally different and conclude that a mea- 
surement of the polarization can be an important source of information to 
uncover the processes behind heavy hadron production. 



2 A+ production in p-p collisions 

In this section we give the prediction of the two component model for the 
cross section as a function of xp- In section 2.1 we show the calculation of 
parton fusion and fragmentation and in section 2.2 the recombination picture 
is presented. 



2.1 A+ production via parton fusion 

In the parton fusion mechanism the A+ is produced via the subprocesses 
QQidO) ~^ cc with the subsequent fragmentation of the c quark. The inclusive 
Xp distribution of the A+ in pp colhsions is given by |TD|, pT| 

daPf 1 ^ f _ , „,.ID 



-y/s J Hab{xa,Xb,Q^) — ^^ dzdp'^dy , (1) 



dxp 
where 

HabiXa,Xb,Q'^) = ^a,b{qaiXa,Q'^)qbiXb,Q'^) 

+ QaiXa, Q'^)qb{Xb, Q^)) -^ \qq 

' dt 

+9a{Xa,Q^)9b{Xb,Q^)-^\gg, (2) 

with Xa and Xb being the parton momentum fractions, q{x, Q^) and g{x, Q^) 
the quark and gluon distributions in the proton, E the energy of the produced 
c-quark and Dj^^/c (z) the fragmentation function measured in e"'"e~ interac- 
tions. In eq. |T], p^ is the squared transverse momentum of the produced 
c-quark, y is the rapidity of the c quark and z = Xp/xc is the momentum 
fraction of the charm quark carried by the A+. The sum in eq. runs over 
a,b = u,u, d, d, s, s. 

We use the LO resuhs for the elementary cross-sections ^ |gq and ^ \gg 

m- 
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where Ay is the rapidity gap between the produced c and c quarks and 
rfi^ = m^ + p^. The Feynman diagrams involved in the calculation of eqs. ^ 
and ^ are shown in fig. |I|. 

In order to be consistent with the LO calculation of the elementary cross 
sections, we use the GRV-LO parton distribution functions |jl2|]. A global 
factor /T ~ 2 — 3 in eq. |l| takes into account NLO contributions [|l^ . 



We take nif. = 1.5 GeV for the c-quark mass and fix the scale of the inter- 



action at Q = 2m^ |T0[. Following M, we use two fragmentation functions 



to describe the charm quark fragmentation; 

D^j,{z) = 5{l - z) (5) 



and the Peterson fragmentation function WA 



N 

Dk./c{z) = — — -T2 (6) 

z[l- l/z-ej{l - z)\ 

with ec = 0.06 and the normahzation defined by J2h I DH/c{z)dz = 1. 



2.2 A^ production via recombination 

The production of leading mesons at low px was described by recombination 
of quarks long time ago |T^. The method introduced by Das and Hwa for 



mesons was extended by Ranft |T^ to describe single particle distributions 



of leading baryons in pp collisions. Recently a more sophisticated version of 
the recombination model using the concept of valons Jl^ has been used to 
study D^ asymmetries in 7r~p interactions |T^. 



In recombination models it is assumed that the outgoing leading hadron 
is produced in the beam fragmentation region through the recombination of 
the maximum number of valence- and the minimun number of sea- quarks. 
Thus A+'s produced in pp collisions are formed by an u and a d valence 
quarks and a c sea quark from the incident proton. Contributions involving 
the recombination of more than one sea flavor are neglected. 
The invariant inclusive xp distribution for leading baryons is given by 

2E da''"' f^p dxi dx2dx3 

p-} =/ F3{Xi,X2,X3)R3{Xi,X2,X3,Xf), (7) 

ay's dXp Jo Xi X2 Xs 

where Xi, i = 1,2, 3, is the momentum fraction of the i^^ quark, F3 (xi, X2, Xs) 
is the three-quark distribution function in the incident hadron and -R3 (xi, X2, 
X3,xf) is the three-quark recombination function. 
The A+ production by recombination is shown in schematic form in fig. 2. 

We use a parametrization containing explicitely the single quark distri- 
butions for the three-quark distribution function 

F3 {Xi, X2, X3) = (3Fu,val {Xl) Fd,val {X2) Fc,sea (a^s) (1 - Xi - ^2 - Xj,)^ (8) 



with Fq (xj) = Xiq (xj) and Fu normalized to one valence u quark. The 
parameters (3 and 7 are constants fixed by the consistency condition 



p L Xi pL Xi Xj 

Fq{Xi) = / dXj dXkF3{xi,X2,X3) , 

Jo Jo 

t,j,k = 1,2,3 (9) 

for the valence quarks of the incoming protons as in ref. ||16||. 



We use the GRV-LO parametrization for the single quark distributions 
in eqs. |] and ^ It must be noted that since the GRV-LO distributions are 
functions of x and Q^, our F3 {xi,X2,X2,) also depends on Q^. 

In contrast with the parton fusion calculation, in which the scale Q^ of 
the interaction is fixed at the vertices of the appropriated Feynman diagrams, 
in recombination there is not clear way to fix the value of the parameter Q^, 
which in this case is not properly a scale parameter and should be used to 
give the content of the recombining quarks in the initial hadron. 

Since the charm content in the proton sea increases rapidly for Q^ grow- 
ing from m'^ to Q^ of the order of some m^'s when it becomes approximately 
constant, we take Q^ = Ami, a conservative value, but sufficiently far from 
the charm threshold in order to avoid a highly depressed charm sea which 
surely does not represent the real charm content of the proton. At this value 
of Q^ we found that the condition of eq. || is fulfilled approximately with 
7 = —0.1 and f3 = 75. We have verified that the recombination cross section 
does not change appreciably at higher values of Q^. 

The value of Q^ in recombination could be different from the scale of the 
interaction in parton fusion. In fact, in the later the scale of the interaction 
should be chosen to be of the order of the hard momentum scale while in 
recombination, the value of Q^ must be chosen in such a way that the true 
content of quarks be present to form the outgoing hadron, allowing in this 
way a more higher value of this parameter. 

For the three-quark recombination function for A+ production we take 
the simple form 



-R3 {xu, Xd, Xc) = a-^^-^-^6 {xu + Xd + Xc- xp) , (10) 

Xp 

with a fixed by the condition Jq dxp{l/a)da'^^^/dxF = 1, where a is the cross 
section for A+'s inclusively produced in pp collisions. Of course, with this 



choice for the parameter a we will not be able to give a prediction for the 
total cross section for A+ production. The recombination model of eq. |^ 
can only give predictions on cross sections relative to the cross section for a 
known process. Thus, we could take a different point of view and chose a in 
order to fit the inclusive xp distribution for a given reaction, e.g. pp -^ p+X, 
and expect that the relative normalizations of the single quark distributions 
account for the desired proportion of another processes cross sections relative 



to the pp ^ p + X cross section, as was made in ref. [jT6[. Instead of it, and 
since in the case of pp -^ A+ + X we don't know the total cross section accu- 
rately, we simply fix a in order to normalize to 1 the A^'s xp distribution as 
given by eq. |^, lossing the possibility to give predictions on the total inclusive 
A^ cross section in recombination. Other forms than eq. [1^ for the recombi- 
nation function have been considered in the hterature [|l7|, |18|, |19|. We have 
observed that the shape of the inclusive cross section of eq. is practically 
insensitive to the form of the recombination function as was pointed out in 
conection with Aq production in pp collisions ||T6[| . 



Using eqs. ^ and |T^ the invariant xp distribution for A+'s produced 
inclusively in pp collisions can be written as 

—= — = alh 2 / dxiFu,vai[xi) 

yjsa dxp Xp Jo 

rxp-xi 
X dX2Fd,val{x2)Fc,sea{xp - Xi - X2) , (U) 

Jo 
where we integrate out over X3. The parameter a is fixed with experimental 
data. 

The inclusive production cross section of the A+ is obtained by adding 
the contribution of recombination (eq. [TT|) to the QCD processes of eq. |T|, 

dxp dxp dxp 

The resulting inclusive A^ production cross section da*°^ /dxp is plotted 
in fig. H a) and b) using the two fragmentation function of eqs. ^ and || 
and compared with experimental data in pp collisions from the ISR 0. As 
we can see, the shape of the experimental data is very well described by our 
model. We use a factor a = 0.92(0.72)/i6 for Peterson (delta) fragmentation 
respectively. 

Although we are not able to explain the normalization of the ISR data 



(we need a global factor of 1000 in eq. |T2]), its xp dependence is well de- 
scribed by our model. 



3 A+ production by IC coalescence and frag- 
mentation 

In a similar approach R. Vogt et al. [^ calculated the A^ production in pp 
and vrp collisions. The two component model they use consists of a parton 
fusion mechanism plus coalescence of the intrinsic charm in the proton. 

In ref. it is assumed that A+ 's can be produced by coalescence of a c 
quark with a ud diquark when the coherence of a \uudcc) Fock state of the 
proton breaks due to the inelastic interaction with the target in pp -^ A+ + X 
reactions. For A+ production in 7r~p collisions, a \udcc) Fock state of the pion 
is considered. 

The frame independent probability distribution of a five particle state of 
the proton is 0, |^ 



by^Ut^^l ...U-O^C 



with N^ normalizing the \uudcc) state probability. 

The intrinsic charm xp distribution is related to P^*" by 
da'"" .. u2 /-i 



dxi 



^ppT^ j dxudxu'dxddxcdxc6 {xp - x^ - x^ - Xc) 

dP^^ 

(14) 



(JjUj iiCajiJU ii' • • • ^"^c 



where the factor /i^/4m^ is the result of the soft interaction needed to break 
the coherence of the Fock state. 0"*^ is the inelastic pp cross section. 

In ref ^ an attempt is made to fix the soft scale parameter at yU^ ~ 
0.2 GeV^ by the assumption that the diffractive fraction of the total cross 
section is the same for charmonium and charmed hadrons. In this way they 
obtain a^^ ~ 0.7 fib in pN interactions at 200 GeV with P^*^ = 0.3%. 

Since we want to compare the prediction of the IC model for A+ produc- 
tion in pp collisions at the ISR energy ^/s = 63 GeV, we prefer to fix the 
unknown product c'^Jlfi'^ with the experimental data. 



The total inclusive cross section for A^ production in pp collisions is then 
given by, 

dxp dxp dxp dxp 

where r represents the fraction of IC production with respect to parton fusion 
(see ref. 0). 

The prediction of the IC two-component model is plotted in fig. ^ c) 
and d) for Peterson and delta fragmentation respectively and compared to 
the prediction of the recombination two-component model and experimen- 
tal data from the ISR 0. We use rhc = 1.8 GeV as quoted in ref. 
and ra^j^fi'^ = 270(320) for delta (Peterson) fragmentation in order to fit 
the experimental data with the IC two-component model. It must be men- 
tioned that, as with the recombination two-component model, the IC two 
component-model can not explain the abnormally high normalization of the 



ISR data and a global factor of 1000 is needed in eq. 15. 



4 The Thomas Precession Model and the A^ 
polarization 

In the TPM the polarization of hadrons is a consequence of the Thomas pre- 
cession during the recombination process ^. 

A A baryon is formed by the recombination of a ud diquark in a spin 
state j = m = and a Q quark originally present in the sea of the projectile, 
so the spin of the A is carried entirely by the Q quark. 

Since the recombining quarks in the projectile must carry a fraction of 
the outgoing hadron's transverse momentum, the Thomas precession appears 
due to the change in the longitudinal momentum of quarks as they pass from 
the projectile to the final hadron. 

In particular, if the A^ is produced through Valence Valence Sea (VVS) 
recombination, since the x-distribution of sea quarks is very step, the A+ 
must get most of its momentum from the valence ud diquark. In this case 
the c sea-quark must be accelerated in passing from the proton sea to the 
produced A+ and the later must be negatively polarized due to the Thomas 
precession. 



On the other hand, if the A^ is produced by the coalescence mechanism, 
since the a;-distribution of the intrinsic charm is hard, the A+ gets most of 
its momentum from the intrinsic c-quark, which must be decelerated in the 
recombination process and consequently the A+ is produced with positive 
polarization ||20[ . 



The amplitude for the production of a A of spin s is proportional to 
{AE + s- oJt), where AE represents the change in energy in going from the 
quarks to the final state in absence of spin effects and cJr~7/(l + 7)-Fx/3 
is the Thomas frequency. 

The polarization in the TPM si given by 

12 [1-3^x^)1 



where 






where mo and iriq are the masses and ptd and pTq the transverse momen- 
tums of the ud diquark and c quark respectively, mA and pta the mass and 
transverse momentum of the Ac, and ^ {xp) = {xq)/xf the average momen- 
tum fraction of the Q quark in the projectile. 

To give a quantitative prediction for the polarization, DeGrand and Mi- 
ettinen take a parametrization for ^ {xp) and obtain a good description of 
experimental data [§]. We explicitly compute ^{xp) [13' 01 using a recom- 
bination model 

(0 (xp) = ^ = ^^"^^^^ffc (18) 

where xp is the Feynman x defined by xp = 2pi^j ^. 

The production of A^ would have a recombination component which con- 
tributes with a smaller fraction to the total cross section. This is a conse- 
quence of the fact that in a proton beam the production of the A^ takes 
place through the recombination of sea quarks only. 

In the production mechanism via intrinsic charm the A^ can only be 
produced by the c-quark fragmentation in a five particle Fock state. If a 
nine particle Fock state is considered uudccuudd) as pointed out in [11 



the Xc distribution would not be as hard as in the five particle Fock state. 
Therefore the recombination would not decelerate the c quark as it happens 

9 



with c quark in the Ac formation via recombination. 

Henceforth, the two mechanisms predict zero polarization for the A^ . 

4.1 The polarization in a two component model 

The polarization is defined as 



Pixp) 



da^ da 



dxp dxi 



I 



da 
dxp 



(19) 



where da^jdxp, dajdxp and da/dxp are the spin up, spin down and total 
cross section respectively. In a scenario where A^s originate in two different 
processes one must take into account the different contributions to the po- 
larization. 

The polarization of baryons in the TPM is the result of the recombination 
process. A+ produced by parton fusion and fragmentation has been shown 
to have a very small polarization pO] which will be neglected here. We have. 



da 



tot 



dal^^ 



da7 



daPf 



We now define the fraction of A+ produced by parton fusion of the total of 
A+s produced as 



hence, from eq. 12 



9{xf) 



daPf 



'daPf 
dxp 


/ 


'da'°'' 
dxp 



g[XF) 



da'''"' 



dxi 



Using eq. ^ we can rewrite eq. ^ 
da 1 



1 — g (xf) dxp 
as 



'rfaf^ 



dxi 



+ 



dal 



g (xp) \ dxp ' dxp J 
and with the help of eqs. ^ and |2^ we obtain the polarization 

dai 



Pixp) = [l-gixp)] 



da\^^ 
dxp 



dxi 



I 



da\^'' dal^^ 



dxi 



or usmg eq. 



16 



P{xp) = -[l-g{xp)] 



12 (l-3e(x^)) 



dxi 



Pta 



(21) 



(22) 



(23) 



(24) 



(25) 



AxM^{l + 3^{xp)f 

where M^ is given by eq. 17. A similar reasoning gives the same eq. 25 for 
the polarization of A+ produced by intrinsic charm coalescence. 
The g (xp) for the two models are shown in fig. ^. 
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4.2 The A+ polarization 

In the coalescence of the intrinsic charm of the proton with quarks in a Fock 



state the xp distribution of the A+ is given by eq. 0. So that we can write 

'' -(..-..) r ^ dx J ^-'v-^--^^' v. (26) 



dxcdxp Jo \Xc + I — Xu — xp 



If the recombination mechanism occurs among the diquark ud in the 
proton and the c quark from the proton's sea then, from eq. |Tl|, we have 

(j^fj rXp—Xc rXp — Xc—x. 



dxcdxp 



pXp—Xc rXp — Xc—Xu 

/ dXu dXdXcF3{Xu,Xd,Xc)R3{Xu,Xd,Xc) . (27) 

Jo Jo 



In order to calculate the polarization of the produced A+'s, we need the mean 
value of ^^c'(^<^'=) (^xp) which is obtained replacing eqs. ^ and ^ into eq. |18 



for each case respectively. The ^{xp) for these two processes are shown in 
fig-i 

We calculate the A^ polarization given by eq. ^ for each production 



mechanism separately using ^ (xp) as given by eq. |T8| with eqs. ^ and |26| 
for recombination and coalescence respectively. 

The A+ polarization is plotted in fig. ^ for the two models; parton fusion 
plus recombination and parton fusion plus IC-coalescence. As in ref. 

we use mp) = | GeV , {Pt)dc = \PTk + (^t) with {k"^) = 0.25 GeV"^ and 
Ax = 5 GeV-^ . We take m^ = 1.5 GeV and m^j = 2.285 GeV . 

5 Conclusions 

We studied the A^ production in pp collisions in the framework of two com- 
ponent models. We used two different versions of the two component model 
and compared them. As we can see in fig. ^ c) and d), both versions can 
describe the shape of the xp distribution for A^'s produced in pp collisions, 
but none of them can describe the abnormally high normalization of the ISR 
data quoted in ref. . At this point it is useful to mention that the normal- 
ization of the ISR data has been object of study in other publications (see 
e.g. ref. ||2^) and none of them was able to explain it. This discrepancy be- 



tween theory and experiment does not exist for charmed meson production, 
which is well described both in shape and normalization with the parton fu- 



sion mechanism plus intrinsic charm coalescence [ 1 1 1 . An analysis of charmed 



meson production using parton fusion and recombination is forthcoming [El . 
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In pp interactions the intrinsic charm component seems to be consistent 



with 1% or less of the total charm sea as stressed in ref. [24]. This suggests 
that the recombination of valence and c-sea quarks may be an alternative 
mechanism to explain the enhancement in the xp distribution at large values. 

The Thomas Precession Model predicts very different results for the A+ 
polarization in pp collisions for the production mechanisms discussed. Unfor- 
tunately, the only known experimental measurement of the A^ polarization 
in pp collisions [§] is not able to give the sign, then we can not obtain defini- 
tive conclusions about the production mechanisms of the A+. However, it 
is interesting to note that the measurement quoted in ref. seems to be 
compatible with a high value for the A^ polarization, despite the large ex- 
perimental errors, which may indicate that the TPM is not adequated to 
describe the polarization of hadrons containing heavy quarks. At this re- 
spect, it is interesting to note that the TPM predicts small values for the 
A+ polarization because of its big mass, (see eqs. |l^ and 0). Moreover, the 
fact that only a fraction of the total number of A^s produced originated in 
a recombination process, supresses the polarization even more. 

Another sensible test that can give information about the production 
mechanisms of hadrons containing heavy quarks is the measurement of hadron 
anti-hadron production asymmetries. In fact, since the parton fusion mech- 
anism predicts a very small asymmetry at NLO, the asymmetry observed 
experimentally in charmed hadron production must be produced by conven- 
tional recombination or the intrinsic charm coalescence mechanism. 

The Ac/Ac asymmetry has been calculated in pp collisions using the two- 
component model with intrinsic charm coalescence [|^. In [0, a factor r = 100 
is needed to describe properly the asymmetry compared to PHYTIA. Such 
very high value for the parameter r seems unnatural because it produces 
an enlargement of the intrinsic charm cross section which can not be easily 
explained. In addition, the intrinsic charm two-component model is not 
able to describe the shape of the experimentally observed asymmetry in D^ 
production in 7r~p interactions ^. 
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Figure Captions 

Fig. 1: Feynman diagrams involved in the calculation of the LO elementary 
cross sections of hard processes in parton fusion. 



Fig. 2: Recombination picture: a c-quark joints a -urf-diquark to build a A 



Fig. 3: a) xp distribution predicted by parton fusion with a Peterson frag- 
mentation function (dashed line) and the sum of the QCD and 
the conventional recombination contributions (solid line). 

b) xp distribution predicted by parton fusion with a delta fragmen- 
tation function (dashed line) and the sum of the QCD and the 
conventional recombination contributions (solid line). 

c) xp distribution predicted by parton fusion plus recombination (full 
line) and parton fusion plus IC coalescence (dashed line) with the 
Peterson fragmentation function. 

d) Same as in (c) for delta function recombination. 

Experimental data (black dots) are taken from ref. 3. 

Fig. 4: Fraction g{xp) of A+ produced by parton fusion and conventional re- 
combination (solid line) and IC recombination (dashed line). In (a) the 
delta and (b) the Peterson fragmentation function was used to describe 
the hadronization. 

Fig. 5: ^{xp) obtained from the conventional recombination model (solid line) 
and from IC coalescence (dashed line). 

Fig. 6: Ajt polarization as a function of xp at px = 0.5 GeV for parton fusion 
plus recombination (lower curves) and parton fusion plus IC coalescence 
(upper curves): 

We show both when a Peterson (solid line) and a delta (dashed line) 
fragmentation function was used. In a similar way for the upper curves 
a Peterson (dotted line) and a delta (dot-dashed line) fragmentation 
function was used in the hadronization. 
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